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Abstract 

In single spin Magnetic Resonance Force Microscopy 
(MRFM), the objective is to detect the presence of an elec- 
tron (or nuclear) spin in a sample volume by measuring 
spin-induced attonewton forces using a micromachined can- 
tilever. In the OSCAR method of single spin MRFM, the 
spins are manipulated by an external rf field to produce 
small periodic deviations in the resonant frequency of the 
cantilever. These deviations can be detected by frequency 
demodulation followed by conventional amplitude or en- 
ergy detection. In this paper, we present an alternative to 
these detection methods, based on optimal detection theory 
and Gibbs sampling. On the basis of simulations, we show 
that our detector outperforms the conventional amplitude 
and energy detectors for realistic MRFM operating condi- 
tions. For example, to achieve a 10% false alarm rate and 
an 80% correct detection rate our detector has an 8 dB SNR 
advantage as compared with the conventional amplitude or 
energy detectors. Furthermore, at these detection rates it 
comes within 4 dB of the omniscient matched-filter lower 
bound. 

1. INTRODUCTION 

Magnetic Resonance Force Microscopy (MRFM) is a re- 
cently developed technique with which physicists can po- 
tentially push the limits of force detection to the single elec- 
tron spin level, with sub-angstrom spatial resolution I1I2I3I . 
The experiment involves detection of perturbations of a thin 
micrometer-scale cantilever whose tip incorporates a sub- 
micron ferromagnet. Any spinning electrons in the sam- 
ple will act as magnetic dipoles, exerting perturbing forces 
that can be measured from cantilever displacements. There 
have been several successful experimental demonstrations 
of MRFM for imaging micron-size ensembles of spins. For 
example, three-dimensional imaging with micrometer spa- 
tial resolution has been achieved |4|. Furthermore, forces 
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as small as 8 x 10 N have been detected using MRFM 
l5l . However, despite several advances, detection of an iso- 
lated single electron spin has not yet been accomplished. 
Progress towards this goal will require advances in physical 
measurements and advances in signal processing of these 
measurements. 

Recently a MRFM method known as OS dilating 
Cantilever-driven Adiabatic Reversals (OSCAR) (6) has 
been proposed to detect single spins. This method, ex- 
plained below, uses a modulated external radio frequency 
(rf) magnetic field to manipulate the electron spins, in or- 
der to produce periodic forces on the cantilever that can 
be detected as small frequency shifts. Detection of these 
frequency shifts identifies the presence of the electron 
spin. If successful, single electron spin detection would 
be an important step towards the long-term goal of three- 
dimensional imaging of subsurface atomic structure \ 7\. 

Unfortunately, accurate single-spin detection in OSCAR 
is hampered by several factors. The spin-induced frequency 
shift signal is extremely weak as a spin induces a frequency 
shift of only one part in 10 4 . Thus long integration times 
are required to detect such a signal. However, spin relax- 
ation and decoherence significantly reduce the usable inte- 
gration time, especially at room temperature. This makes 
the use of cryogenics (cooling the experimental apparatus 
down to a fraction of a degree Kelvin) necessary to reduce 
these effects. However, in this regime the measurements 
are prone to thermal noise from various sources. Noise and 
decoherence effects must be taken into account by the de- 
tection algorithm to achieve the most accurate and reliable 
single spin detection. Very simple detectors are the base- 
band "amplitude detector" and "energy detector" which op- 
erate on a frequency demodulated version of the cantilever 
position signal. Such detection schemes are widely used in 
MRFM, NMR spectroscopy, MRI, and other applications. 

In this paper, we present a new approach to baseband 
detection in OSCAR experiments. The detector is based 
on a random telegraph model for the baseband signal in- 
corporating Poisson-distributed random spin reversals, ran- 
dom initial spin polarity, and Additive White Gaussian 
Noise (AWGN). To accurately decide between the hypothe- 



ses of spin absence and presence, we propose a hybrid de- 
tection scheme which combines optimal Bayes and Gen- 
eral Likelihood Ratio (GLR) detection principles imple- 
mented with Gibbs sampling. Simulations show that our 
proposed detector can significantly outperform the conven- 
tional baseband amplitude and energy detector for realistic 
post-demodulation signal-to-noise ratios (SNR). 

The outline of the paper is as follows. After briefly re- 
viewing the OSCAR experiment in Sec. 2, we describe the 
proposed signal detection scheme in Sec. 3, and present re- 
sults of numerical simulations in Sec. 4. 
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2. DESCRIPTION OF EXPERIMENT 



Fig. 1. Schematic of the MRFM experiment. 



Fig- HI is a schematic description of the OSCAR exper- 
iment. In OSCAR, a submicron ferromagnet is placed at 
the tip of a cantilever which sits at a distance of approx- 
imately 50 nanometers above a sample. In the presence of 
an applied rf field, electrons in the sample undergo magnetic 
resonance if the rf field frequency matches the Larmor fre- 
quency. Since the Larmor frequency is proportional to the 
field from the magnetic tip, and because the tip field falls 
off rapidly with distance, only those spins that are within a 
thin "resonant slice" just the right distance from the tip will 
satisfy the condition for magnetic resonance and contribute 
to the signal. 

If the cantilever is forced into mechanical oscillation by 
positive feedback, the tip oscillation induces small shifts 
in the Larmor frequencies of the spins. Specifically, the 
tip motion gives rise to an oscillating magnetic field which 
sweeps the Larmor frequency of the spins in the resonance 
slice back and forth through resonance. This causes the spin 
to reverse polarity synchronously with the cantilever mo- 
tion, and in return, the spin reversals affect the cantilever 
motion by changing the effective stiffness of the cantilever. 
When an electron spin is present the spin-cantilever inter- 
action can be detected by measuring small shifts in the pe- 
riod of cantilever oscillation using laser interferometric can- 
tilever position sensing. For more details about OSCAR, 
see 171 151191 . 

A classical (non-quantum) electro-mechanical descrip- 
tion of the spin-cantilever interactions can be developed in 
the framework adopted by Berman et al | 10 1 and Rugar et 
al 1 1 1 1. We briefly review this framework here. Consider a 
spin in a rotating frame which rotates at the frequency of the 
applied rf magnetic field, Bi (Fig.|2}. The effective mag- 
netic field B eJ y (t) in this frame is given by 



B eff (t) = Bxi + &B {t)k, 



(1) 



where i and fc are unit vectors in the x and z directions in 
the rotating frame, B\ is the amplitude of the rf magnetic 
field, B a (t) is the magnetic field caused by the magnetic 
tip of the cantilever, and AB a (t) = B Q (t) — (J r f/l is the 



off -resonance field magnitude (7 is the gyromagnetic ratio). 
Spins are in resonance on the spherical shell (resonant slice) 
defined by those spatial locations for which ui r f matches the 
Larmor frequency jB (t). 

If AB a (t) varies sufficiently slowly such that the adia- 
batic criterion 



dAB (t) 
Jt 



« iBl 



(2) 



is met, the spin can be assumed to remain aligned or anti- 
aligned with the vector B eJ y (t). This is the spin-lock condi- 
tion. Let the vertical position of the cantilever tip be denoted 
by z where z — denotes its rest position. Under the influ- 
ence of the external rf field Bi (t), electron-spin forces, and 
random (thermal) force noise F n (t), the motion of the can- 
tilever tip obeys the simple harmonic oscillator equation: 



mz(t) + Tz(t) + kz(t) = 



\»\G 2 z(t) 



y/G*z(t)* + B\ 



(3) 



where m is the cantilever's effective mass, k is the can- 
tilever spring constant, T is the friction coefficient char- 
acterizing cantilever energy dissipation, \fi\ is the magni- 
tude of the spin magnetic moment, G = dB oz /dz is the 
z-direction field gradient at the spin location. The natural 
mechanical resonance frequency of the cantilever is given 
by uo — y/k/m, and T can be related to the cantilever 
quality factor, Q, via V = k/ {lo Q). 

Under the small tip displacement approximation \ Gz\ <C 
B\ , we obtain 



mz(t) + Tz(t) + (k + Ak)z{t) 



F n (t), (4) 



where Afc = — \/j,\G 2 /Bi. This shift in spring constant 
results in a shift Au of the cantilever resonant frequency: 
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In a version of the protocol called "Interrupted OSCAR," 
the Bi field is turned off every T s ki P seconds over a 
half cycle duration (ir/ui ) to cause periodic transitions 
between the spin-lock and anti-spin-lock spin states (see 
middle panel of Fig. |4}. In the spin-lock state the spin 
aligns with the field B eJ y(i) and in the anti-lock state the 
spin aligns with — B ej gr(t). Therefore, the frequency shift 
Aui of the cantilever alternates between the two values 
±iw (|/x|)G 2 / (kBi) with period T s ki p . In the absence of 
noise (F n (t) = in (|3}) the cantilever motion can be ex- 
pressed as the frequency modulated (FM) signal: 



z(t) = Acos[uj t+ \ s(t')dt' + 



(6) 



Here A is the cantilever oscillation amplitude, 8 is a random 
phase, and s is equal to if no spin coupling occurs, while 
it is equal to a periodic square wave of period 2T s ki P and 
of amplitude \Aui \ if spin coupling occurs. Thus, in this 
ideal noiseless case, the presence of spin coupling can be 
perfectly detected either by detecting a spectral peak near 
|Aa> | radians in the periodogram or by frequency demod- 
ulation of z to baseband (incorporating subtraction of the 
known center frequency lo ) followed by amplitude detec- 
tion, energy detection, or other algorithm, as discussed be- 
low. As baseband and narrowband are equivalent represen- 
tations we focus on the baseband method here. These meth- 
ods correlate the baseband signal against the known square 
wave signal derived from Bi. The resultant signal, which 
we call y(t), forms the statistic which is used for spin de- 
tection, as illustrated in Fig. [5] 
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Fig. 2. In the coordinate system rotating at w r /> the off-resonance 
field AB , and therefore the effective field B e g it), vary with time. 
Under the spin lock assumption, the spin always follows ±B ej gr(£). 



Unfortunately, in a practical (non-ideal) experiment the 
cantilever tip-displacement signal z(t) is degraded by sev- 
eral factors which reduce correlation peak detection accu- 
racy. One factor is the presence of laser interferometric 
measurement noise. This adds a noise floor to the demod- 
ulated square wave signal s(t). Another factor is spin re- 
laxation which over a period of time causes the spin to 
go out of alignment with the effective field B ej y. While 
several models for single-spin relaxation have been pro- 
posed I1 11101 . a full understanding of the physics of single- 
spin relaxation interactions with cantilevers remains open. 
One model, adopted here, is that the single spins main- 
tain spin lock or anti-lock states but spontaneously and asy- 
chronously change polarity during the course of measure- 
ment at some rate A flips/second. In the sequel we develop 
an optimal single-spin detection approach under a random 
Poisson model for these polarity flips. 

3. SIGNAL DETECTION IN NOISE 

The signal detectors we will consider operate on the 
baseband output signal y(t) of the frequency demodulator, 
e.g. a Phase-Lock Loop (PLL), followed by multiplication 
by a square wave reference pit") G {±1} of period 2T s ki P , 
whose transitions are synchronous with the (known) rf turn- 
off times (see Fig. [3}. 
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Fig. 3. Baseband detector frequency demodulates the interfer- 
ometric signal, correlates the output against a square wave p(t) 
whose transitions are synchronous with the turn-off times of the 
rf field Bi(t), and generates a test statistic, e.g. accumulated 
squared frequency deviations, for detecting presence of a spin. 



We model the baseband output y(t) of the frequency de- 
modulator and correlator as a random telegraph plus addi- 
tive Gaussian white noise (see lower panel of Fig. [4}- Let 
[0, T] be the total measurement time period over which the 
correlator integrates the measurements, and let {ri},i = 
1..N, be the time instants within this period at which spin 
reversals occur. We assume {Ti] are the arrival times of 
a Poisson process with intensity A. Consequently N is a 
Poisson random variable with rate AT 1 12 1. Thus, y(t) = 
s(t) + v(t) where v(t) is AWGN with variance = er 2 , and 



s(t) is a random telegraph signal containing only the ran- 
dom transitions: 



s(t) 



\Aoj 



N 

£( 

i=0 



-!)*<?(- 



f - 



-l - n 



(7) 



where <j> is a random variable that takes on ±1 with equal 
probability, representing a random initial spin polarity, tq = 
0, tn + i = T, and g(t) is the standard rectangle function: 
g(t) — 1 for t £ [0, 1] and g(t) = otherwise. If there are 
no random spin flips in the time period [0, T], then s(t) = 
0|Ao; o | is constant over time, which we obtain in by 
using the convention that when N = 0, t d = and Ti = oo. 

The baseband spin detection problem is to design a test 
between the two hypotheses: 



H (spin absent): y(t) = v(t) 
Hi (spin present): y(t) = s(t) ■ 

fort £ [0,T]. 
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Fig. 4. Top: Sample cantilever position signal, z(t), at 10 kHz. 
Middle: Sample rf magnetic field magnitude, B\, has synchronous 
half-cycle skips at 1 ms, 2 ms, and 3 ms. Bottom: In the presence 
of a single spin, s in Eq. \(frhas both deterministic transitions due 
to the rf skips at 1 ms, 2 ms and 3 ms, and random ones due to 
spin relaxation. The random transitions, {Ti}, occur as a Poisson 
process. The initial polarity is (j> — 1 for this example. The noisy 
signal at bottom is s with AWGN contamination. 



Conditioned on the random parameters {t^}, N, <j>, the 
signal s(t) is deterministic and known. Under this condi- 
tioning the optimal detection structure would be the simple 
matched filter 1131 



y(t f )s(t',cl>,T,N)dt' 
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< 



(9) 



where s(t;<f>,T,N) is a synthesized random telegraph 
signal of the form parametrized by |Aw | (assumed 



known), (f>, r and N. The value 77 is a threshold that can ei- 
ther be set to satisfy a probability of false alarm {Pf) con- 
traint Pp < a, a £ [0, 1], or as a function of the prior 
probabilities a \tx[P{Hq) / P{H\ )]+b where a, b are known 
constants. In the former case the detector is called the most 
powerful (MP) test of level a, which has maximum proba- 
bility of detection (Pd), while in the latter case the detector 
is called the minPe detector as it achieves minimum average 
probability of decision error (minPe). 

As the values of the random parameters are always un- 
known, we call the detector (|9j the omniscient matched fil- 
ter, which is unimplementable. However, as the omniscient 
matched filter is optimal for known parameter values it es- 
tablishes a useful upper bound on performance. 

Perhaps the simplest baseband detection scheme, and the 
most widespread in MRFM applications, is the amplitude 
detector which acts as if there were no random flips and 
declares a spin present if the magnitude of the average am- 
plitude of the correlator output exceeds a threshold 
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(10) 



where r] is a threshold set to give the desired Pp. Improved 
performance can be obtained by explicitly accounting for 
the equally likely initial polarity and assuming AWGN to 
derive the minPe detector. The amplitude detector d 1 Oft is 
the minPe detector under the assumption that y(t) is a ran- 
dom polarity constant embedded in AWGN. This is a valid 
assumption when there are no random spin flips over the 
integration period [0, T). 

When there are random spin flips over [0, T] due to spin 
relaxation and decoherence, the performance of the ampli- 
tude detector suffers. Indeed, as the number of random flips 
increases the average amplitude of y(t) converges to zero. 
As the energy of s(t) is independent of number of transi- 
tions, transition times, and initial polarity, it is natural to 
propose an energy detector lfl3l 
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where 77 is a threshold set to give the desired Pp. It can 
be shown that the energy detector is a minPe test for the 
case that v(t) is additive white gaussian noise, s(t) = 
Aw cos{2-Kt/T s ki P + 9), and 6 is uniformly distributed over 
[0, 2n] [ 14 1. It can also be shown that the energy detector is 
the minPe test under a Gaussian approximation to the ran- 
dom telegraph process in the limit of high SNR 1 15 1. Note 
that in our case the detection performance of the energy test 
is independent of the average flip rate A. 

As we will show in the sequel, the performance of the 
amplitude and energy detectors can be far from the optimal 
performance achieved by the omniscient matched filter de- 
tector. 



3.1. The Hybrid Bayes/GLR Detector 

The minPe detector for a signal with random parame- 
ters is a Bayes likelihood ratio test that averages an omni- 
scient likelihood ratio test statistic over all random parame- 
ters fl4l : 



logA(y) 

e t ,n [^[/(yjT.Jvr,^!)]] 



(12) 



= log- 



f(y\H ) 



Hi 
> 
< 
Hn 



As above r\ is a threshold selected to achieve a de- 
sired level a of Pp. The function / is the joint p.d.f 
of {y(t)}te[o.T] parametrized by the random parameters 
t, N, <p, and B x [-|yl] denotes conditional expectation over 
random variables x given event A. 

While the expectation over 4> in d!2t is simple to eval- 
uate, the expectation over {{t^}, N} is very difficult since 
the integration region is of very high (infinite) dimension. 
An alternative to performing this second expectation is to 
invoke the Generalized Likelihood Ratio (GLR) principle. 
The GLR consists of replacing the unknown parameters in 
(I12> by Maximum Likelihood (ML) estimates. 
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where, again, 77 is a threshold chosen for a desired Pp. Note 
that in Jl 31 we have averaged over <p while we have max- 
imized over {{ri},N}, leading to what we call a hybrid 
Bayes/GLR test. 

As y(t) is a conditionally Gaussian random process 
given {t,;} and N, the log-likelihood function in d!3> can be 
simplified by invoking the Cameron-Martin formula 1 16 1: 



is guaranteed to converge to the maximizing solution, is 
to more efficiently search over the space by Gibbs Sam- 
pling II 811191 . As we know the Poisson intensity A, we can 
generate samples {{t^}, N} from the prior Poisson distri- 
bution so as to maximize the log-likelihood function. As 
these samples are more likely (on the average) to mimic the 
actual behavior of the parameters, we obtain a reduction in 
search complexity. 

The general description of the Gibbs sampler is as fol- 
lows. Supposed there is a random vector variable X = 
[xi, X2, ■ ■ ■ , x p ] T having density function /x from which 
we want to sample. Suppose also that we can simulate the 
i-th element of X given samples (already simulated) of the 
other elements: 

X% \x\ j X2 1 • ■ • 1 ^i — l j Xi-\-\ , . . . , Xp 

~ fi(Xi\x!,x 2 , ■ ■ .,Xi-i,x i+ i, ...,x p ) fori = l..p 

(15) 

Then a Markov sequence, x'*' = [x± , . . . , Xp] T , can be 
simulated by the recursion 

Xt 1] ~ h( Xl \4\...,x^), 
Xt 1] ~ h{x 2 \x^\ x f,...,xV), 



fp{ x p\%l \ x 2 ; * * • j 1 )■ 



(16) 

After a certain amount of burn-in time X^*^ 1 , t > Tf,, will 
have stationary distribution /x ■ In our case, since the arrival 
times {r;} are generated from a Poisson process, the condi- 
tional distributions fl!6l > are easy to sample from, because 
they are conditionally uniform. 

4. SIMULATION METHODS AND RESULTS 



logA(y) = max < log cosh 

T,N 
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y{t)a + (t;r,N)dt 
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(s+(t;r,N)) 2 dt 



(14) 



where s + (t;T, N) is the synthesized telegraph wave Q 
having initial polarity cf) — 1 and parametrized by r and 
N. It is well known that for a sufficiently large integration 
time T the minPe and GLR tests are identical (see for ex- 
ample 1 17 1). Thus we can assert that the hybrid Bayes/GLR 
test Jl 3i is an asymptotically optimal test. 

3.2. Solution via Gibbs Sampling 

The maximization in dl4> by exhaustive search over the 
uncountably infinite dimensional space of possible parame- 
ter values, {{t^}, N}, is impractical. An alternative, which 



The objective of our first three simulations was to com- 
pare the detection performance of the matched filter, the am- 
plitude detector, the energy detector, and the Bayes/GLR 
detector on the basis of Receiver Operating Characteristic 
(ROC) curves, which are obtained by empirically generat- 
ing the pairs (Pp , Pp> ) for each detector. In our simulations, 
the four decision rules l|9}, i ll Oi l. ( II It . and Jl 3I > were used 
to generate the ROC curves in the Matlab 6.1 environment. 
Based on the Monte Carlo methodology 1 18 1, we generated 
samples {y^\n)}, yd(n) = y(nT s ), under both Hypothe- 
sis and 1, where T s was the sampling period. The samples 
were input to the detector being evaluated, and Pp, and Pp 
were statistically calculated. 500 detection trials were per- 
formed under each hypothesis. For each ROC curve, the 
above process was repeated with a range of decision thresh- 
old values 77. This range of thresholds was chosen to ade- 
quately sample the domain Pp e [0, 1]. 



The simulation parameter values were chosen according 
to typical OSCAR experimental values. The signal dura- 
tion T was 3 s and the sampling period T s was 0.5 ms. The 
rest of the parameters were set to: k — 1 x 10~ 3 N/m, 
w /(2tt) = 1 x 10 4 Hz, B x = 0.2 mT, G = 2 x 10 6 
T/m, and /i = 9.3 x 10" 24 J/T. With these parameters 
the (noiseless) signal amplitude, |s(f)|, was 0.928 Hz ac- 
cording to 0. Two values of A, the average number of 
random flips per second, were evaluated. The detector 
noise was assumed AWGN and the noise variance was ad- 
justed to investigate the effect of SNR, which is defined as 
Wlog 10 [(l/T) \s(t)\ 2 dt/a 2 v }. 
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Fig. 5. Simulated Receiver Operating Characteristic (ROC) 
curves for the matched filter, energy detector, amplitude detector, 
and hybrid Bayes/GLR detector, at SNR = -25 dB and A = 1 event- 
per-second. Unlike the other detectors, the matched filter assumes 
complete information on the parameter values and is not imple- 
mentable. 



We ran the Gibbs sampler for 5,000 iterations for the hy- 
brid Bayes/GLR detector. Fig.|3]and[6]show ROC curves for 
SNR = -25 dB and -20 dB, respectively, for A = 1 event- 
per-second. In both cases, our hybrid Bayes/GLR detector 
significantly outperformed all the other detectors except for 
the unimplementable matched filter. The matched filter had 
complete information about the random flip times, and as a 
result it achieved almost perfect detection for both SNR val- 
ues. In Fig.0the value of A was increased to 10 events-per- 
second and the SNR was held at -20 dB. As expected, the 
performance of the amplitude and hybrid Bayes/GLR detec- 
tors degrades, while the matched filter and energy detector, 
whose performance does not depend on A, are not affected. 
In Fig. U the power curves for all detectors are plotted as 
a function of SNR for A — 1. Here all detectors perform 
at the same false alarm rate Pp = 0.1, and we can make a 
quantitative SNR comparison by fixing the detection perfor- 
mance level at Pp> = 0.8, say. To attain this detection level, 
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Fig. 6. Simulated Receiver Operating Characteristic (ROC) 
cur\'es for the matched filter, energy detector, amplitude detector, 
and hybrid Bayes/GLR detector, at SNR = -20 dB and A = 1 event- 
per-second. 



the energy detector and amplitude detector require SNR's of 
at least -14 dB and -17.5 dB, respectively, while the hybrid 
Bayes/GLR detector only requires -26 dB. As compared to 
the amplitude detector, this represents an improvement of 
almost 9 dB in SNR performance using our proposed detec- 
tor. Furthermore, the performance of our proposed detector 
is only 4 dB worse than the performance bound of -30 dB 
established by the matched filter for this level of Pp and 
Pry. Note that the amplitude detector outperforms the en- 
ergy detector for low SNR's but not for high SNR's. This is 
explained by the fact that even though the energy detector is 
not affected by random flips, at low SNR's its test statistic 
is dominated by the noise. 

In another simulation, we investigated the role of the 
number of Gibbs samples on performance of the hybrid 
Bayes/GLR detector, shown in Fig. [9] It is evident that 
performance improves as we increase the number of Gibbs 
samples. For example, at Pp = 0.1, Pp> increases from ap- 
proximately 0.35 to 0.65 if we increase the number of Gibbs 
samples from 100 to 500. It increases further to around 0.9 
and 0.95 if 500 or 5,000 Gibbs samples are used, respec- 
tively. Such improvements in performance are significant 
but yield diminishing returns as the number of Gibbs sam- 
ples is increased beyond 500. 

5. CONCLUSION 

In this paper we presented a hybrid Bayes/GLR approach 
to detecting the presence of single spins for the OSCAR 
MRFM experiment. We have shown by simulation that the 
Bayes/GLR detector performs significantly better than the 
classical amplitude and energy detectors. The improvement 
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Fig. 7. Simulated Receiver Operating Characteristic (ROC) 
curves for the matched filter, energy detector, amplitude detector, 
and hybrid Bayes/GLR detector, at SNR = -20 dB and A = 10 
events-per-second. 
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Fig. 8. The power curves ( Pd vs. SNR ) for the four detectors 
studied in this paper for Pf =0.1 and A = 1 event-per-second. 
At Pd =0.8 the hybrid Bayes/GLR detector performs within 4 dB 
of the bound established by the matched filter. 



in detection performance is due to the fact that, unlike the 
classical detectors, the new detector estimates the unknown 
values of the random spin reversal times and the initial po- 
larity. Of course this performance improvement comes at 
the price increased implementation complexity. This com- 
plexity increases in the random reversal rate A due to the 
necessity to perform Gibbs sampling over an increasingly 
large number of probable spin reversal sequences. Nonethe- 
less, for the experiments in the first simulation with A = 1 
event-per-second, the run time of our detector (5,000 Gibbs 
samples) was only on the order of about half a minute per 
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Fig. 9. Juxtaposition of ROC curves of hybrid Bayes/GLR detec- 
tor, obtained with different numbers of Gibbs samples in the max- 
imization step, at SNR = -20 dB and A = 10 events-per-second. 
Performance improves as the number of Gibbs samples increases. 



3-second measurement record (our code was implemented 
in Matlab 6.5, under WindowsXP on a 2.26GHz PC, with 
510MB RAM). 

An interesting extension of our results would be to as- 
sume that the frequency shift |Acj | is also unknown. This 
would lead to a hybrid Bayes/GLR detector which detects 
the peak over the spectrum of the signal in addition to max- 
imizing over the number and positions of the transitions. 

The hybrid Bayes/GLR detector was derived using a 
baseband signal model consisting of a random telegraph 
wave with Poisson transitions and AWGN. This signal 
model is theoretically justified under the spin-lock assump- 
tion. The validity of the spin-lock assumption remains to 
be established. More sophisticated signal models of the 
cantilever measurements, and associated detection methods 
which bypass frequency demodulation and operate directly 
on those measurements, are currently under investigation. 
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